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Abstract 

A novel approach is presented in this paper to analyze attitude precision and 
control for an instrument gimbaled to a spacecraft subject to an internal disturbance 
caused by a moving component inside the instrument. Nonlinear differential equations of 
motion for some sample cases are derived and solved analytically to gain insight into the 
influence of the disturbance on the attitude pointing error. A simple control law is 
developed to eliminate the instrument pointing error caused by the internal disturbance. 
Several cases are presented to demonstrate and verify the concept presented in this paper. 

Keywords: Attitude dynamics and control, feedback control, spacecraft dynamics 

Introduction 

This paper was originally motivated by developing a fine-pointing control 
algorithm for enabling the GIFTS (Geostationary Imaging Fourier Transform 
Spectrometer) instrument to meet its mission requirements [1], GIFTS combines a 
number of advanced technologies to observe atmospheric weather and chemistry 
variables in four dimensions. It will enable meteorological soundings equivalent to those 
achieved by simultaneously launching of 16,384 closely spaced radiosonde balloons 
within a 600-km diameter circle about every 10 seconds. It will revolutionize 
atmospheric science and meteorological forecasting. GIFTS in s trument/mi s sion affords 
an opportunity to space-qualify a significant number of new technologies, payload and 
non-payload specific, for future generation remote sensors. 
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The GIFTS attitude precision for the instrument pointing requires the use of a 
lightweight, low power, inexpensive, and highly reliable mini-star-tracker for estimation 
of the spacecraft orientation. Note that star cameras are among the most attractive attitude 
sensors, because they provide three-axis attitude information with high accuracy. 
Considerable work has been done in developing the mini-star- tracker, along with attitude 
estimation algorithms producing significant achievements as documented in Refs. [2-10], 
The main role of the mini-star-tracker that will be integrated into an Optical Pointing and 
Stabilization Control System (OPSC) is to achieve the GIFTS pointing requirements. 

The goal of the OPSC is to orient and stabilize the GIFTS optical instrument with 
respect to the Earth. OPSC will scan the earth as well as the near-stellar region off the 
earth’s horizon in a step-stare procedure. Imaging occurs during a ten second stare at a 
specific point on earth, followed by a step motion of one FOV to a new target within one 
second, covering the disc of the earth approximately every 15 minutes. OPSC will also 
stabilize the optical boresight against disturbance motions generated by the spacecraft 
that carries the GIFTS optical instrument. 

The original strategy of the GIFTS OPSC is shown in Figure 1 . The optical 
platform is mounted on a two-axis flex-pivot gimbal for boresight pointing control while 
vibration is reduced through a high bandwidth jitter control mirror. Attitude 
measurements are made by a dual star camera system and jitter is measured by an MEMS 
rate sensor integrated with the star cameras. Gimbal and jitter mirror angle measurements 
are made using high precision electro-optic sensors. A separate processor is used for 
system control. 

This paper begins with a brief introduction of two-body dynamics. All 
differential equations of motion derived for many differential cases are based on the 
fundamental concept of two-body dynamics. A simple configuration of a rigid body 
with an internal moving part representing the GIFTS instrument rigidly attached to a 
spacecraft will be discussed. Parametric studies will be performed to understand how the 
moving part influences the overall system motion, in particular the inertial pointing. The 
other configuration to be studied consists of one rigid body (spacecraft) gimbaled by 
another rigid body (instrument) with a moving part inside. Three rigid bodies are 
involved in this configuration making the case extremely complex in deriving equations 


2 



of motion for designing a control-torque law to meet the instrument attitude precision 
requirements. Furthermore, it is very difficult, if not impossible, in performing 
parametric studies for general cases. For simplicity, without losing generality, we will 
limit the parametric and numerical studies in one rotational axis. 


| Line of Sight 
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Figure 1: Optical Pointing and Stabilization Control System (OPSC). 

Basic Formulation 

The basic configuration consists of the spacecraft and the optical instrument with 
a moving mirror. OPSC is designed to orient and stabilize the optical instrument by 
providing a control torque/force against the spacecraft. During the control maneuvering, 
coupled motion of the spacecraft and the instrument takes place. To describe the overall 
system motion, multiple-body dynamics are involved. Basic formulations for two-body 
dynamics will be briefly described in this section [11], 

From Figure 2, let m a be the mass of the body a (for our case, the spacecraft that 
carries the optical instrument), \ a be the velocity of the coordinate origin 0 ; , (o a be the 

angular velocity of the body, and p a be the distance vector from 0 to an arbitrary point 


3 


in the body. Subscript a signifies the associated quantity for the body. The linear 
momentum of the body over the domain Q a is described by 


P a = f Q (v a + (o a x p a ) dm a = m a (v a + co a x c fl ) 


( 1 ) 


where 


c 


a 



I P , dm 
jQ„ a 


is the definition of the center of mass (CM). 



Figure 2: Body and inertia coordinates for two-body dynamics. 

Similarly, define mt as the mass of the body b, \ h as the velocity of the 
coordinate origin fl ft , to 6 as the angular velocity of the body b, r p as the distance vector 
from 0 to fl A where the optical instrument is maneuvered for fine pointing, and p, as the 
distance vector from 0 to an arbitrary point in the body b. The linear momentum of the 
body b over the domain Q h is 
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( 2 ) 


P b = L K+tOjXp h ]dm b 

= L [ v a + ® a x ( r „ + Pi) + *o p X P J dm b 

= m h [v„ + to a x (r, + c 4 ) + (o p x c 6 ] 
where the angular velocity ( 0 ; of body b relative to body a is 

(Op=(O b -(» a 


and the center of mass c h is 


c h = — L p/, dm b 

m h Ja b 


The basic momentum equation of motion about point 0 (a or b ) is 

h = | Q P x ( v » + v ) dm (3) 

with v = p being the velocity of an arbitrary point in Q . Equation (3) is applicable for 
both bodies a and b. For body a, replace h, p , and v 0 by h a , p a , and v fl . For body b, 
replace h, p , and v o by h A , p b , and \ h . The angular momentum of the body a about 
its reference point 0 a is 

h , = Pa X (■ v a + to , x p )dm a = m a (c a xv) + 1 ,to fl (4) 

where 


l a«>a = | Q P , X (CD„ x p a )dm a = | (p>„ 1 - p ,p' )dm a 


to 


and 1 is a 3 by 3 identity matrix. Note that the second equality is written in matrix form 
for easy numerical implementation. Similarly, The angular momentum of the body b 
about its reference point 0 is 

K = L P/,x(v, + to,xp b )dm b 

= L Pb X [v a + to„ X (r + p 6 ) + © x p b ]dm b (5) 

= m h {c h xv a ) + l ba (o a +l h (o p 
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where 


!/,»,, = j n Pft xOo,, X P/ )dm b = J Q (p[p 6 1 - p/ p' )dm b 


to, 


and 


I 6«®« = L P* x t®, x ( r P + P b Wm b 

j p [Pi (* p + P*)l “ (Fp + P,)p[ ]<K 


0) , 


= [I b +m b (c T b r p l-r p c T b )](o a 
The temporal derivative of h from Eq. (3) becomes 

h = \ Q P x (v 0 + y)dm + \ a p x (v. + y)dm 

= li2 VX ( V ° +V )^ m + | i2 P X ^f (6) 

= pxv 0 + r 

From Eqs. (1) and (2) , the equations of motion for translation are: 

Pa=f.-f, (7) 

and 

P,=f,+f p (8) 

where f, and f/, are external forces applied to a and 6 respectively, and f p is the internal 
force acting at the joint 0 ,. From Eq. (6), the equation of motion for rotation about the 
origin 0 for the body a is: 

h „+T x P„=T-^+r ? xf p (9) 

and the equation of motion for rotation about the joint 0 ; , for the body b is 

ha+K+c^xrpjxpj =T b +r p (10) 

where T a and t b are external toques applied to a and b respectively, and i p is the 

internal torque acting at 0 ;/ 

Summing Eqs. (7) and (8) for system translation yields 

P a +Pi = f „+ f z, => P= f (11) 


The rate of change in the total system momentum equals to the total external force. On 
the other hand, adding Eq. (9) to Eq. (10) for system rotation gives 
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h + v a xp = T 


( 12 ) 


where the total angular momentum 


h = K + h * + r p x P/, = m (c x v a ) + I(0„ + 1 


,C0 

ab p 


and the torque is 


T = T + T, + r xf, 

a b p b 


The total moment of inertia over £2 and £2, about 0 is 

a b a 


1=1 + L, K r /> + P*) 7 ( t p + Pa) 1 - ( r „ + P/,)( r p + P b Y¥m b 
= I a +I b + m b (rj r p 1 - r/ p ) + m b (2r T p e h 1 - r/ b - e b r T p ) 


(13) 


The center of mass for the whole system over £2 u and Q b is 

me = m B c„ + m b (c 4 + r ;7 ) ; m = m a + m b 
The mix moment of inertia is 

Kb = [P l (r p + P/,) 1 - P h ( r P + P bY ¥m b =l b + m b (rjc, 1 - c 6 rj ) 

Equation (12) implies that the rate of change in the total system angular momentum 
equals the total external torque. The equation of motion to orient and stabilize body b 
can be obtained by inserting Eq. (2) into Eq. (10) to yield 

K = (v.+^xr^xKcjXK+wJ + Tj+T^ (14) 

Equations (11), (12), and (14) are differential equations of motion for a two-body 
dynamical problem with one body to be reoriented for fine pointing. For easy numerical 
implementation, they may be reformulated in terms of matrix form as 

p = -K,x]p-f 


h = T-[(o„x]h-[v a x]p (15) 

K = {[hiXl + fTJv* +[w«x]r,)x]Kc J x]}(T A +® p ) + t 4 + T, 

where each vector quantity has three components, and T ha is the transformation matrix 
from the body frame a to the body frame b. For any vector o> , its [ox] is 



CO, 


0 

- CO 3 

co 2 

0) = - 

co 2 

>=> [tox] = 

CO 3 

0 

-CO, 


CO, 


-co 2 

CO\ 

0 _ 
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Based on Eqs. (1), (2), (5), and (13), the momentum equation in matrix form is 



P 


ml 

-[wzcx] 

-T UK c a x] 



£? = M v => 

h 

= 

[?«cx] 

I 

^ ab 


(16) 


_v 


[m b c b x]T ba 

ha 

h 

_°V 



where 


1 r 

c = — 
m 


m a c a+m b (T ab c b +r p ) 


; m = m + m. 


1 = 1, + W 'lb + ™ b [(rjr, )1 - r/ p ] + m b [2 (rj T 6 c A ) 1 - v/ h Y ab - T 6 c,r 


^ab ~ T \ b | 


X b +m b 


*~ba 


'Kvd^vX]} 


r. = {0"i[Kv,)*-VXlt 


'-ba 


Note thatp and h contain components in the body frame a, but h ; contains components in 
the body frame b. Angular velocity becomes 

v = jvr 1 £? (17) 

To compute the transformation matrix from one frame to the other involves 
kinematics. Let the quaternion vector q be the four components expressed by 


01 

02 

03 

04 


(18) 


The quaternion vector equation is 


and its scalar equation is 


d l . 

— q = —(q to-wxq) 
dt 2 * 4 


d 1 
— q = — 0 ) q 
dt 4 2 M 


(19) 


( 20 ) 


Equations (19) and (20) in matrix form are 


0i 


02 

1 

03 

~2 

_04 _ 



0 

to. 

-to. 

CO, 

01 

-0 3 

0 

(0 ] 

co 2 

02 

ry 2 

-0) s 

0 

CO, 

03 

-a»i 

-to. 

-CO, 

0 

_04 


( 21 ) 



Note that both bodies have the same form of kinematics. For the spacecraft, replace to 
by to a and q by q a . For the instrument, replace to by to A and q by q 6 . The 
transformation matrix from inertia frame to body frame becomes 


l BI 


2 2 2 2 

q,-q 2 -q, +q 4 

2 (q 2 q l -q 4 q 2 ) 
2 (g i q i + q 4 q 2 ) 


2 (q t q 2 +q 4 q 4 ) 

2 2 2 2 

-q; + q:~q:+q; 

2 (q 3 q 2 -q 4 q,) 


2 (g l q i -q 4 q 2 ) 

2 (q 2 q 2 +q 4 q 2 ) 

2 2 2 2 

-q;-q 2 +q 2 +q; 


( 22 ) 


This matrix and its transpose are used to transform vectors from the inertia frame to the 
body frame, and vice versa. 


One-Rotation-Axis Plane Motion 

For simplicity, let us assume that the moving mirror is the only moving part of the 
spacecraft together with the optical instrument. In other words, the optical instrument is 
locked to the spacecraft as shown in Figure 3. 



Figure 3: Body and inertia coordinates for one -rotation-axis plane motion. 


The spacecraft linear momentum with the total mass m a and the coordinate origin 
at the center of mass is 


P a = m a v a =m a (Xe x + Ye Y ) 

= m a ( X cos 6+ Y sin 6 > )e x + m a (-X sin 0+ Y cos 0)e y 


(23) 
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where 6 is the rotation angle from the inertia frame and the body frame. The unit vector 
e x indicates the direction of the inertia axis X, and e Y is the unit vector for the inertia 

axis Y. The quantities X and Y are the speeds of the coordinate origin along e x and e Y , 
respectively. The instrument linear momentum with the mass m b is 


P, =»h 


v a + o) a x (r, + c b )] = m b [(Xe x + Te Y ) + 6e z x t e y ] 


= m b [(X cos 6 + Y sin 6-1 6)e x + (-X sin 6 + Y cos 6)e ] 


(24) 


where the unit vectors e x , e y , and e z give the directions of the body coordinates x,y, and 

z, respectively. The quantity 6 is the angular velocity of the spacecraft, and t is the 
distance from the center of mass of the spacecraft to the center of mass of the instrument. 
The moving-mirror linear momentum (assume point mass) is 

Pv = rn 3 (v fl + <0 a X r s + s) = m s (Xe x + Te Y ) + 6e z x (fe y + .ve x ) + ve x 


= m s [(; 


(25) 


(X cos 6 + Tsin 6 - 1 6 + s)e x + (-X sin 6 + Y cos 6 + s6)e y 

where s is the distance from the center of mass of the instrument to the moving-mirror 
that is considered as a point mass.. 

Note the following coordinate transformation 

x= X cos 6+ Y sinf? 

[y = -X sin 6+Y cos 6 
Equations (23), (24), and (26) reduce to 

P„ = ™ a [xe x +ye y ] 

P b =m b [(x-£6)e x +ye y ] (27) 


(26) 


p s = m s [(x-£0 + i)e x + (y + s6)e y ] 

Spacecraft angular momentum about the reference point o is 

h " = L, p “ X (v “ + 0) “ X 9 • ) dm “ = I " t0 “ = (28) 

Instalment angular momentum about the reference point o is 
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( 29 ) 


h„ =ro i c 4 xv a +I ia a> a 

= m/e y x[xe x +ye y ] + l ba 0e z 
= [-m/x+\ ho 0]e. 


Moving-mirror angular momentum about the reference point o is 
h = m> [rx(v„+o)„xr v + r)] 

= m s [ i e y + se x ] x |[x - i 6 + i ]e x + [y + s(9]e y | 
= m s {-£[x- id + s] + s[y + sd]}e z 


Translational equation of motion for the whole system is 


P=P a +Pb+Ps =° 


| (m a + m h +m s )Jc- (m b + m s ) i 6 = -m s s 
(m a + m h + m. )j + m s s6 = 0 


Rotational equation of motion for the whole system with p = 0 is 

h=h a +h b +h v = 0 

=» ~{m b + m s ) tx + m s sy + [I fl + l ba + m s ( t 2 + s 2 )]6 = m s is 


Equations (31) and (32) yield the following matrix equation of motion 


m a +m b + m s 

0 

~(m b +m s y 

X 


-m s s 

0 

m a +m b + m s 

m s s 

y 

= 

0 

_ ~(m b + m s )l 

m s s 

l a+ l ba+ m s(t 2+s2 )_ 

0 


m/s 


Define the non-dimensional quantities as 

_ x _ y _ s 
x = y =— , s =- 


//i| — , rn^ — 

m + m , + m ~ m + m, + m 

a b s a b s 

j l a+ l ba+m / 2 

{m a +m h + m s )f 

then the matrix equation of motion, i.e., Eq. (33), becomes 


1 

0 

-m t 

X 


-m 2 s 

0 

1 

nys 

y 

- 

0 

-m ] 

m 2 s 

Ij + m 2 s 2 

0 


m 2 s 


(30) 


(31) 


(32) 


(33) 


(34) 


(35) 


Equation (35) may be solved to yield the following differential equations to be integrated 
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(36) 


x = 


-m 2 [ij - m 1 + (1 - m 2 )m 2 s 2 1 s 
I, - m[ + (1 - m 2 )m 2 s 2 


y 


9= 


-(1 -m\)m 2 s s 


I, - m\ + (1 - m 2 )m 2 s ‘ 
-(1 - m, )m 2 s 


Ij — ml + (1 -m 2 )m 2 s 2 
Analytical solutions for Eqs. (36), (37), and (38) are 

(l-7?7,)tan _l s 1 - m 2 ) m 2 j (l, - m [ ) 


9 = 


y = 


^{\ x -ml){\-m 2 )lm 2 

(l - 777 j ) m 2 log [ij - 777 2 + (l - 777 2 ) 777,5 2 1 


2 ( 1 - 777 ,) 


and 


(37) 

(38) 


(39) 


(40) 


X = 777 ! 9 — 777,5 


(41) 


Equation (39) shows the relationship among the rotation angle 9 , the moving-mirror 
motion J , and the inertia ratio Ij with given mass ratios 777, and 777, . The angle error 9 
induced by the moving-mirror motion is proportional to J for sufficiently small ~s and 


large 1, , because under such condition tan 


-1 


5 V ( 1 _ ' 77 2 ) 777 2 /( I 1 - 777 1 2 ) 


approaches 


s ^(l - m 2 )m 2 /(l l - ml j . Figure 4 illustrates the complex but interesting relationship 

where the moving mass, the spacecraft mass, and the instrument mass were given from an 
earlier GIFTS design. For large moving-mirror motion ( 3T > 1 ) and relatively small 
inertia ratio ( I, < 20 ), the pointing-angle error may be larger than 0.02 degree. The 
inequality Ij < 20 implies that the moving mirror is located away from the center of mass 
of the spacecraft ( i » 0 ) and the moment of inertia ( I fl + \ ha + m s l 2 ) is relatively small 
such that the ratio of l a + \ ba + mj’ 2 to [m a + m b + m s ) i 2 is less than 20. The relationship 

surface is quite nonlinear in the upper-front-end comer. To minimize the pointing-angle 
error via a control torque requires a nonlinear control law that will be developed in the 
following section. 
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Moving-mirror mass m s = 10 Kg 
Spacecraft mass m a = 726 Kg 
Instrument mass mb = 190 Kg 


Angle 

error 

(deg) 



Normalized moving- 
mirror motion s 


Inertia 
ratio Ii 


Figure 4: Relationship among the angle error 6 , the moving-mirror motion s , and the 
inertia ratio I, . 


Two-Rotation-Axis Plane Motion 


Let the instrument be allowed to rotate along the z-axis. Let co a = 0 a be the 


rotational speed of the spacecraft about the z-axis and co b 


0 b be the rotational speed of 


the instrument also about the z-axis. Figure 5 illustrates the quantities required for the 
following development of the dynamic equations. 

First note the coordinate transfonnation from the body coordinates of the 
spacecraft to the inertial coordinates 


e x 

e Y 


cos 6 a 

- sin 9 a 

V 

sin0 a 

cos 6 a 

1 

a j 

i 


(42) 


and the transformation from the body coordinates of the instrument to the inertial 
coordinates 



cos 6 b 

- sin Q b 

e x’ 

sin G h 

cos 6 b 

_ e y’_ 


( 43 ) 
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Figure 5: Body and inertia coordinates for two-rotation-axis plane motion. 


The spacecraft linear momentum is 

p =m(\ +» xc ) 

tr a a \ a a a / 


= 


( X - c a cos 0 a e a )e x + (Y-C a sin 6 a 6 a )e Y 


(44) 


where \ a is the velocity at the coordinate origin o (i.e., the contact point of the spacecraft 
and the instrument), X and Y are its components along e x and e Y in the inertial frame, 
respectively. For simplicity, without losing generality, the spacecraft center of mass c a 


is assumed to be in the direction of e. . The instrument linear momentum is 


P/, =m i [v a +o i xcJ 


= m. 


(* 


)e x + m b ( 


X-c h cos G h G b e x + m,, ( Y - c k sin 0, 0,. | e 


b^b ) Y 


(45) 


Note that the direction of the body coordinate e y , is chosen to pass through the 

instrument center of mass. The linear momentum of the moving-mirror traveling in the 
direction e x , is 


14 


P, =w,[v a +o 4 x(c 4 + s) + s] 

= m s \^X - (e, cos 6 h + s sin G h )G b + cos G h s e x 
+ m s Y + ( s cos 6 h - c b sin Q b )G b + sin G b s e Y 


Spacecraft angular momentum about the reference point o is 


K = ip., P “ X[V “ +0> “ X P^ dm a = M o ( C a X V a ) + I a °>« 

= \_- m „ C a ( X C0S 6 a + Y Sin 6 a ) + T A C Z 

Instalment angular momentum about the reference point o is 

h* = f P, x [ v fl + to, x p, ] dm b = m h (c, x y fl ) + 1,(0, 

= \- m b c b ( x cos 0 b + Y sin O b ) + r , G b e z 

Moving-miaor angular momentum about the reference point o is 
h , = m s ( c b + s) : X [v a + CO, x (c, + s) + s] 

= m s j-(c, cos O h + s s i n ®b) x ~ ( c b s i n @b ~ s cos + ( c l + s 2 )4> ~ c b ^\ ( 


To perform the parametric study, let us define the following non-dimensional 


quantities 


X — —■ Y — — • v — — • r = — • r = — • 

x*- ? 1 ? ° ? '-'a ? L 'h i 


m h + m s 

/«, = — ; m 2 = 

m a + m b + m s 


m a +m b +m s 


[ m a +m b +m s 


y ; 2= i 


m,+m b +mAc 


p h +h = h +l b h 

[m a +m b + m s )c (m a + m b +m s )c 


+ m b + m s I c 


where c is an arbitrary positive constant intuitively set to be c = | c a | + c, | . Normalizing 
the total linear momentum p = p + p, + p s yields 

Px X - (1 - m l )c a cos Q a Q a - (m,c, cos G b + m 2 s sin G b )G b + m 2 s cos O h 

P- = . ‘ “ . (52) 

P y \ Y - (1 - m, )c a sin O a O a - (m,c, sin G h - m 2 s cos G b ) G b + m 2 s sin G h 
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Normalizing the sum of the angular momentums h/, and h s produces 

h + h v = |-[m 1 c 6 cos 9 h + m 2 J sin 9 b ]X - [m t c h sin 0 h - m 2 J cos d h j\Y 
+ [I 2 + m 2 c; + m 2 J 2 ]6 h - m 2 c b s }e z 

Furthennore, normalizing the total angular momentum of the system h = h„+ h/,+ h v 
yields 


( 53 ) 


h = |-[(1 - m { )c a cos 0 a + mf b cos 6 h + m 2 s sin d h \X 

- [(1 - m l )c a sin 9 a + mf h sin 9 h - m 2 J cos 9 b )]Y (54) 

+ l A + t 1 2 +m i c l +m 2 s 2 ]9 h -m 2 c b s}e z 

The angular momentum in Eq. (53) or (54) has only one component along the Z axis. 


Differential Equations of Motion 

With no external force applied to the system and zero initial condition, differential 
equations of motion for the overall system translation are 


P = 


P x 

Py 


0 

0 


=>p = 


P x 

Py 



P x (t = 0 ) 
Py(t = 0) 


0 

0 

o' 

0 


(55) 


where p is the non-dimensional quantity defined in Eq. (50) with the assumption of zero 


initial linear momentum, i.e., P(f = 0) = 0 . The differential equation for the overall 
system rotation with no external force and torque is 

h = T-v fl xp = 0 => h = 0 (56) 
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The differential equation of motion for the instrument rotation is 

h / J +h v=-V a X(p Z) +p s ) + T p 

= - x { 0,, X [ (m b + m s )c h + m s s] + m s s} + r f 


= | r p - m s Xs sin 0 b + m Ys cos 0 b 

+ X0 b [~m s s cos 0 h + (m h + m s )c b sin 0 b ] 
-Y6 h [m s s sin 9 h + (m h + m s )c b cos <9 6 ]} e z 

that produces the non-dimensional equation of motion 


h, +h v = r p - m 2 X s sin <) h + m 2 Y s cos 0 b 


+ X0 b -m 2 s cos 0 h + ni\C b sin 0 h 


-Y0 b m 2 s sin 0 h + m x c b cos 0 b e 


where 


( m a 


+ m b + m 


V ' ( 


m, + m b + m ) c 




r D = 


( 57 ) 


(58) 


i m a 


+ m b + m I c 


.)< 


(59) 


In view of Eqs. (52), (53), and (54), the whole system can be described by the 
following four equations [see Eq. (15)] 

p x = 0; p } = °; K = o; 

ti_ = Y p - m 2 X J sin 0 h + m 2 Y J cos 0 b + X 0 h [~m 2 J cos G b + m^c b sin Q b ] (60) 


-Y 9 b [m 2 J sin 9 b + mf b cos 9 h ] 

where h_ and fi_ are the components of h and h, +h s in the inertial direction e z . 
Equation (60) can be integrated to solve for the time histories of p x , p, , h, , and ti_ . 

Note that the third equation in Eq. (60) is valid only when the condition p x = p y = 0 is 

satisfied at all times and no external torque is applied. At any time t, the quantities X , Y , 
9 a , and 9 b can be updated by using Eqs. (52), (53), and (54) to first form 

£? = Mv + p (61) 
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where 



for use in computing a new quantity fi_ via integrating Eq. (60). 

On the other hand, one may prefer to have a conventional set of equations of 
motion that involve physical accelerations. Differentiating Eqs. (52), (53), and (54) and 
then applying Eqs. (55), (56), and (58) produce the matrix equation of motion 

Mx = f (66) 

where M is defined in Eq. (63), and 

1 " 

Y 

(67) 

0 a 
0 b _ 
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and 


-m 2 s c 0 b + c a (m i - \)0 2 a s 0 a + 2 m 2 s 0 b s 0 b - 0 b (c b m t s0 b - m 2 s c0 b ) 


f = 


-m 2 s s 0 b - c a (m l - \)0l c 0 a - 2 m 2 s 0 h c 0 h + 0 b (c b m l c0 b + m 2 s s 0 h ) 

c b m 2 s +m 2 sX0 b c0 b +c a (m x -l)X0 a s0 a -c a (m x -\)Y0 a c0 a 
-c b m x X0 h s0 b +c b m 1 Y0 b c0 b -m 2 7lYc0 b +2s0 b -Xs^] + Ym 2 J0 b s0 h 


( 68 ) 


T p -2 m 2 s s0 h +c b m 2 s 


What is the torque r p in Eq. (68) that will drive the instrument to satisfy the precision 
pointing requirement? This question will be answered below. 


Control Torque for Attitude Precision Pointing 

The equality 0 h = 0 h = 0 represents the condition where the instrument has no 
angular velocity relative to the inertial frame. The differential equations of motion for 
the instrument rotation with 0 h = 0 h = 0 are 


h b + h s = -m x c b X + m 2 sY - m 2 c b s 
h + h, = r + m^Y 


(69) 


Differentiating the top equation of Eq. (69) and substituting it into the bottom equation 
yields the toque required for the instrument fine pointing 


z p = -m x c b X + m 2 sY - m 2 c b s 


(70) 


where s is a pre-specified quantity; X and Y may be measured by using accelerometers 
Equation (70) shows that the control torque is a weighted sum of the acceleration J for 

the moving mirror and the accelerations X and Y of the spacecraft at the joint with the 
instrument attached. Note that other control techniques, such as predictive control [12], 
may also be used for fine-pointing the instrument subject to known and/or unknown 
periodic disturbances. 
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For the case where the spacecraft is sufficiently large in comparison with the 
instrument such that 

X = 0 and 7 = 0 (71) 

the torque shown in Eq. (70) will then approach to 

r„ = -m 2 c b s (72) 

System translation for 9 h = 9 h = 0 is derived as follows. Substituting 0 h = 0 h = 0 
into Eq. (52) with the aid of Eq. (55) yields 

X - (1 - m l )c a cos 9 a 9 a + m 2 s - 0 (73 ) 

7 - (1 - m, )c a sin 9 a 9 a = 0 (74) 

From Eqs. (54) and (56), the total angular momentum for 0 h = 0 h = 0 is 


-[(1 - m l )c a cos 0 a + nif b ]X -[(1 - m l )c a sin 9 a - m 2 s ]Y + \9 a - m 2 c b s = 0 


Equations (73), (74) and (75) produce the following differential equations to be 
integrated 


~ m i ( J i - c a (i - m l ) { c b cos 0 a + sin 9 a [ -m 2 s + (1 - m, )c a sin 9 a ] } ) 5 
- (1 - m l ) 2 c] + (1 - m x )c a (m 2 s sin 0 a - m x c b cos 9 a ) 


i 2 s sin 9 a - m x c b cos 9 a 


e = (\-m x )m 2 {c h - c a cos 0 a ).y 

Ij - (1 - m, ) 2 cl + (1 - m x )c a ( m 2 s sin 9 a - m x c b cos 9 u ) 

Integrating Eqs. (76), (77), and (78) simultaneously and analytically is quite difficult, if 
not impossible. A numerical example is given in the following section to illustrate the 
concept developed in this paper. In particular, we will show how the control torque 
works for the instalment precision pointing and its influence to other quantities such as 
the system translational and rotational displacements. 
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Parametric and Numerical Analyses 

A representative case will be studied in this section. Given certain parameters, 
such as the mass and inertia ratios, we will derive the control torque for instalment 
precision pointing and show its relationship with the mirror motion and its influence to 
the overall system motion. 


Control Torque with Given Mass and Inertia Ratios 

Let us assume that the mass-center offset ratios of the spacecraft and the 
instrument away from the joint o of the two bodies are given by. 

c=- = - 0.9; cj= — = 0.1; c = |c fl | + |c A | (79) 

c c 

where c is the total mass-center offset. In addition, assume that the mass ratios m, and 
m 2 among the spacecraft mass m a , the instrument mass m b , and the moving-mirror mass 
m s , and the inertia ratio I, related to the spacecraft inertia are given as 


m b + m s 


m a +m b + m s 


= 0.1 ; m 2 = 


m. 


+ m,. + m. 


= 0.01 


I, = 


I 


= i; i 2 = 


L 


(80) 


= 0.1 


(m a +m b + m s )c 2 ’ 2 (m a + m b +m s )c 2 

These ratios imply that the instalment mass plus the moving-miaor mass is 1 0% of the 
total system mass, whereas the moving-miaor mass is only 1% of the total mass. The 
spacecraft inertia about the point o is equal to the inertia contributed by the total system 
mass and the total mass-center offset. 

Differentiating Eqs. (76), (77), and (78) to solve for X , Y , and d a , and 
substituting X , Y , and Q a into the resultant equation yield the control torque required 


for instrument precision pointing 

_ _ [2.767O5 + O.328O5(Tsin20„-cos20„)]Txl(T 4 
Tp ~ -0.3439 + 0.0081(5 sin0 fl -cos 0 a ) 

[2.25633(0. 1 + 0.9 cos 6 a ) 2 (s cos d a + sin 6 a )] 5 2 x 1 0 7 
[-0.3439 + 0.008 l(j sin d a - cos 0„)] 3 


(81) 
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Given a mirror-motion profile such as the square wave that yields s , s , and s , and the 
spacecraft angular displacement 0 a , the control torque time history can be computed 
from Eq. (81). In practice, one would use Eq. (70) with the insertion of measured 
quantities X and Y , and pre-specified s and s to compute the torque r . With the 

control torque computed and applied to the system, the translational and rotational 
displacements are described by Eqs. (76), (77), and (78). For parametric analysis, let us 
rewrite Eqs. (76), (77), and (78) in terms of s as the independent variable, rather than 
the time t, to yield 

dX ~ m 2 ( J i - c a (1 - m, ) { c b cos 0 a + sin 0 a [~m 2 s + (1 - m, )c a sin 0 a ]j) 

— =" — = = — ; = ( 8 ^ 

ds I; - (1 - m, )~c a + (1 - m 1 )c a (m 2 s sin 0 a - m^c b cos 0 a ) 


dY (1 f m 2 c a (c^- c a cos 9 a ) sin 0 a _ 

ds Ij — (1 — m l ) 2 c a 2 + (1 - m l )c a ( m 2 s sin 0 a - m { c b cos 0 a ) 


(83) 


ddg Jl-ffl,)m 2 (c, -c^cosOj 

ds I; - (1 - m l ) 2 cl + (1 - m l )c a (m 2 s sin 0 a - m x c b cos 0 a ) 


(84) 


Integrating Eqs. (82), (83), and (84) simultaneously for 1 > s > 0 produces the following 
three relationships, i.e., X versus s shown in Figure 6, y versus s shown in Figure 7, 
and 0 versus s shown in Figure 8. 



Figure 6: System displacement on X-axis versus moving-mirror displacement 
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Figure 7: System displacement on Y-axis versus moving-mirror displacement 


x 1 0 3 



Figure 8: Angular displacement on Z-axis versus moving-mirror displacement 


Observe that Figure 6 and Figure 8 show a linear relationship for small angle 9 a . 

Assume that the mirror is moving as a sine wave with amplitude a and frequency 
co such that 

s - a sin cot 

s = aco cos cot = aco\J\ - sin 2 cot = acOyjl-(s /a)' (85) 

s = - aco 2 sinryt = -co 2 s 

Substituting Eq. (85) into Eq. (81) and plotting the result f p against s and 0 a in the 
range of 0.1 > s > 0 and nj 4 > () a > 0 with co = 1 Hz yields Figure 9, showing that large 
amplitude J and displacement 9 a will require higher torque T p for fine-pointing control. 
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Figure 9: Control torque r versus moving-mirror displacement s and spacecraft 
angular displacement 0 a 


Control Torque for 0 a - 0 

Consider another case where 9 a = 0 . System equations (76), (77), and (78) for 
linear and angular velocities become 


X 

l 

-m 2 [\ -c a c b (\-m,)]s 

y 

0 


Y 


K 


(1 -mfm 2 (c b -c a )J 


Y = !i - 0 - m i f c a ~ 0 - n h K c b m i 


( 86 ) 


Differentiating Eqs. (76), (77), and (78) to solve for X , Y , and 0 a with 0 a = 0 for linear 
and angular accelerations yields 


X 

Y 


m 


[\~ c a c b(\- m i)] s c] ( i - m ] f ml (c a - c h ) 2 s .v ‘ 


r y 

c a {\-m x f m l(c b -c a ) 2 J 2 

f 


(1 - m, )m 2 (c 6 - c a ) j c a (1 - m, ) 3 /n 2 3 (c a -c t ) 2 si 2 


r 




(87) 
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The control torque for inertia pointing when 9 a = 0 becomes 


T - -mf b X + m 2 s Y - m 2 c b s 

_ c b [ J i ~ gfl ~ m S] C 1 ~ 

Ii — (1 — m i f c] ~ 0 - K c h m x 
| Cg (g fl - c b f [i, ~ c] (1 - Wi ) 2 ] (1 - n\ ) 3 miss 2 
[i, - (1 - m l ) 2 c 2 a - (1 - m l )c a c b m x ]‘ 


( 88 ) 


Equation (88) implies that the control torque r p can be expressed in terms of s , s , and 


s at any given angle 9 a , because the linear accelerations X and Y are purely induced 
by the mirror motion with the absence of external forces. Inserting the numerical values 
shown in Eqs. (79) and (80) into Eqs. (86), (87), and (88) yields 

X -3.07102 xl0~ 2 j 

Y = 0 (89) 

9 a [ 2.55682x10 ~ 2 s _ 

-3.07102xl0" 2 5 -1.2185x10 ^jj 2 

-5.29523 x 1 O' 4 s 2 (90) 

2.55682 x 10“ 2 5 + 1 .50433 x 10‘ 4 n 2 

and 

r /; =-6.92898 xl0" 4 7-5. 17388 xl0“ 6 5 5 2 (91) 

It is clear that the acceleration J dominates the dependent quantities to be computed in 
Eqs. (90) and (91). 



Numerical Simulation 

Assume that the moving-mirror is traveling as a sine wave such that 

J = a(l - cos cot); J = aco sin cot; J = aco 2 cos cot (92) 

Let the amplitude a and the traveling frequency co be 

a = 0.1; co = 1Hz (93) 
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Integrating Eq. (66) with the insertion of numerical values defined in Eqs. (79) and (80) 
and zero initial condition generates the time histories for X , Y , 0 a , and 0 h shown in 

Figure 10 and for X , Y , Q a , and 0 b shown in Figure 11. The control-torque time 

histories are given in Figure 12. The curves marked in red in Figure 10 through Figure 
12 show the time histories with no control torque. The curves in green give the time 
histories with partial control feedback only from the acceleration ~s of the mirror motion 
[see Eq. (72)]. The curves in blue represent the time histories with full-control feedback 

as defined in Eq. (70) from the acceleration J , and the accelerations X and Y of the 
spacecraft at the joint with the instrument attached. 

Intuitively, we may consider only the feedback of the original acceleration source, 
i.e., J due to the mirror motion, and ignore other accelerations induced by the source. 

The time history in green for the pointing error 0 h shown in Figure 12 illustrates a clear 

reduction with control half that of the time history in red without control. Nevertheless, it 
is seen from the lower part of Figure 12 that the control torque, marked in green, with the 

absence of X and Y feedback expresses an overshoot in magnitude relative to the full- 
control torque in blue. It reflects the displacement overshoot, i.e., from positive 
displacement in red without control (cross the blue line on the zero axis) to become 
negative displacement in green with mirror acceleration control only. With full control 

including J , X , and Y feedback, the pointing error 0 b marked in blue vanishes 
completely. 

The time histories for the system displacement X with and without control are 
approximately two orders in magnitude larger than those for the displacement Y . The 
same statement is also true for the velocity time histories. This is due to the fact that the 
simulation begins with the angle at 0 a = 0 and the mirror motion initially perpendicular 

to the Y axis. In addition, since the instrument including the moving mirror is much 
smaller in weight and inertia than the spacecraft, the control torque required for the 
instrument pointing does not induce much into the spacecraft rotation angle 0 u and its 
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angular velocity 6 a . Asa result, the control torque would influence the X-axis motion 
much more than the F-axis motion. 




0 a 

(deg) 


0.3 

0.2 

0.1 

0 


Time (sec) 


Time (sec) 



0 b 

(deg) 


Figure 10: Time histories for X, Y , 6 a , and () h with and without control torque 




Time (sec) 




Y 


0 h 

(deg/ sec) 


Figure 1 1 : Time histories for X , Y , Q a , and 0 h with and without control torque 


27 



Time (sec) 

Figure 12: Time histories for pointing error 0 h and control torque z p . 

Concluding Remarks 

Parametric and numerical studies are performed in this paper to analyze attitude 
precision dynamics and control for an instrument gimbaled to a spacecraft. We focus on 
the attitude pointing error caused by the internal disturbance from a moving mirror inside 
the instrument. First, we examine the relationship among the pointing angle error, the 
inertia ratio, and the moving-mirror displacement, assuming that the instrument is locked 
to the spacecraft. An analytical solution is derived for the complex, but very useful, 
relationship that provides a road map for an engineer to determine how to design an 
optimal configuration for an instrument with a moving mirror. Secondly, we examine the 
case when active control is needed to counter balance or reject the internal disturbance 
for attitude precision control. A simple control law is developed producing the control 
torque that is capable of eliminating the pointing error induced by the internal disturbance. 
The control torque is a weighted sum of the mirror moving acceleration and the system 
linear acceleration at the point where the instrument is attached to the spacecraft. A 
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system/control engineer may minimize the control torque by reducing the weighting 
coefficients that mainly depend on the instrument configuration. Note that all analytical 
solutions in this paper are derived under the assumption of single rotation axis. The 
concept and approach are believed to be applicable to general cases. Numerical results 
from a sample case indicate the importance of an appropriate control torque otherwise 
overshoot may occur that, in turn, waste the control energy. 
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